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ABSTRACT 

The low coupling efficiency of semiconductor laser emissions into single mode 
optical fibers place a severe restriction on their use. Associated with these 
conventional optical coupling techniques are stringent alignment sensi- 
tivities. Using holographic elements, the coupling efficiency may be 
increased and the alignment sensitivity greatly reduced. Both conventional 
and computer methods used in the generation of the holographic couplers are 
described and diagrammed. The reconstruction geometries used are shown to be 
somewhat restrictive but substantially less rigid than their conventional 
optical counter-parts. Single and double hologram techniques are examined 
concerning their respective ease of fabrication and relative merits. 
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HOLOGRAMS FOR LASER DIODE - SINGLE MODE 
OPTICAL FIBER COUPLING 


I. INTRODUCTION 


The use of optical fibers for the propagation of light emitted from a 
semiconductor laser (or laser diode, LD) has been an area of intense study for 
the past 5-10 years. Single mode LDs, in which only the fundamental mode is 
emitted, are proving to be a most attractive source for high data rate optical 
communication systems. When the output from a single mode LD is coupled into 
a single mode optical fiber, fiber dispersion (the spreading of the pulse 
shape due to differing propagation speeds for the modes in the fiber) is 
eliminated. Another benefit in using single mode optical fibers is that the 
emerging beam is well-behaved spatially, having a Airy disk far-fleld 
diffraction pattern. 

In the past five years, a great deal of effort has been expended in 
decreasing the fiber transmission losses. With these losses now on the order 
of 0.2 dB per kilometer, emphasis has been turned to the large losses 
associated with coupling the LD output into the fiber (1-4). 

The techniques currently employed for the coupling of LD light into a 
fiber consist of the following: 

1. Butt-end joint: The fiber is cleaved so that its front face is 

' perpendicular to the fiber core. It is then butted 

up against the LD. 
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2. Tapered fiber; Some of the fiber’s cladding is etched off, by 

chemical or mechanical means, exposing the fiber’s 
core. It is then positioned near the front facet on 
the LD. 

3. Fiber-tip microlens: The tip of the fiber is melted an4 formed into a 

hemispherical lens. This end is then positioned 
such that the LD is at the lens’ focal point. 

4. Conventional optics: Lenses (typically 2 crossed cylindrical lenses) 

are placed in front of the LD to conform its 
output beam and focus it onto the fiber’s core. 

These techniques provide coupling efficiencies of less than 40%^. (The 
40% figure has been reported in technical journal articles a laboratory 
results. A reproducable coupling efficiency of 20% is considered quite good 
(5)). The causes of this poor coupling efficiency are primarily attributable 
to the small acceptance angle of the fiber's 5 micron core (Numerical Aperture 
of about 0.2), and the large non-symmetrlc divergence of the LD beam (20 and 
80 degrees parallel and perpendicular to the 0.7 micron by 10 micron LD lasing 
region. All of these techniques also result in very stringent alignment 
tolerances for the positioning of the fiber and any other optical elements, 
with respect to the LD front facet. Alignments of approximately 1-2 microns 
perpendicular and parallel to the LD front facet must be maintained for losses 

*The coupling efficiency is given by the ratio of the amount of power 
launched into the fiber and the LD's output power. 
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of lesi? than 1 dB. The spacing between the LD and fiber directly Impacts 
coupling. The alignment tolerances In this case are In the 20 micron range 
for <1 dB loss. 


II. WAVEFRONT CONVERSION COUPLER 


A different approach to this problem Is to examine the wavefronts of LD 
and fiber In the Fraunhofer diffraction region. At some plane In between the 
two devices these fields can be superimposed resulting In an Interference 
pattern. Theoretically, if some sort of wavefront converting structure could 
be constructed and placed In this plane to effectively match the two 
"antennas" (the LD and the fiber), a coupling efficiency of 100% could 
result. This procedure would actually convert the wavefront diverging from 
the LD Into the Airy disk far-fleld pattern of the fiber. This type of 
wavefront conversion falls Into the general category of holography. So If a 
hologram representing the wavefront converting structure would be placed In 
this plane and then Illuminated by the LD emission used to construct the 
hologram, all of the LD light Impinging on the hologram would have Its 
wavefront changed Into the fiber’s Fraunhofer diffraction pattern. The 
converted wavefront would then be totally accepted by the fiber core. 

In considering diffraction from an aperture, the Fraunhofer region Is 
defined as the distance, Z, from the aperture where the phase change over the 
aperture is assumed to be much less than 1 radian. If the aperture Is 
considered to be rectangular with the maximum value of Xj^^ + yi^ * (d/2)^ 
where d is the linear dimension of the maximum extent of the aperture, then 
this phase restriction implies: 
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( 1 ) 


or z » /X 


( 2 ) 


Single-mode optical fibers have core diameters of 4-8 microns. LDs' 
maximum front-facet linear dimension Is typically 10 microns. For the fiber 
and LD, the Fraunhofer assumption Implies, assuming light with wavelength 830 
nanometers . 


"f 



77 microns 


z 


LD 



120 microns 


Thus, If the hologram’s positioning plans; Is chosen to be 150 microns 
away from the LD and fiber, the Fraunhofer diffraction condition Is valid. 


The amplitude distribution in the Fraunhofer diffraction pattern of the 
circular fiber core is the Airy disk pattern given by: 


U(r„) 


exp (jkz) exp 




(3) 


where 1 is the radius of the fiber core; r Is the radius coordinate 


In the 


observation plane; Z is the distance from the aperture (fiber core) to the 
observation plane; and Jj^ Is the first order Bessel function. 
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Similarly* the electric field a distance Z away from the LD's front facet 
is given by the Fraunhofer diffraction pattern for a rectangular aperture.* 




exp (JKx) exp M 


2 2 ' 
2z 


1 1 
X y 


»(*o.y„) .=n> 



(5) 


where 1^^ is the linear dimension of the aperture in the x direction; ly is the 
linear dimension in the y direction. 


With the elements positioned as in Figure 1, the electric field at the 
plane P is given by the superposition of the two incident fields and U^. 


■ + 


u, 


The resultant intensity distribution is; Ip ■ 


U„ 


Up* I 


*Equations (3) through (5) were derived assuming plane waves incident on the 
aperture. The field incident on the front facet of the single-mode LD (TEj^q 
mode) is roughly planar, as in the case for a rectangular waveguide (6). 
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For the case of a monochromatic field: 


Ip 


FAiSe IS 

OiF POOR quality 


-lu.nC*,. y„>l^+ 2 Mx,,y„)a(x„.-/„)cos(nx„.y„) 


'LD' o' 'o‘ 
“ 0(^0 . yo» 


Investigation of the terms yields; 



( 7 ) 


( 8 ) 



( 9 ) 


2 2 

X + y 
o ''o 


Z 2 

X + y^ 
o “^o 


'f((Xo.yo)”^(%’yo^ " — Sz — ^ 

( 10 ) 
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Figure 2 shows a cross section of Ip (xo, yo ■ HX z/21y) with X ■ 830 nra, z 
150 microns, ly"0.7 microns, lx*10 microns, 1“6 microns* These plots were 
generated by the computer program PF.FOR (listed In Figure 13). 

III. COUPLING EFFICIENCY 

The coupling efficiency, n can be decomposed Into two principal 
terms; the transfer efficiency n j., and the intrinsic coupling efficiency 

n 1« 


n c " ’I t 1 

depicted in Figure 3. 

The transfer efficiency involves the amount of input radiation absorbed 

by the hologram and Its diffraction efficiency. It has been shown that for 

thick holograms transfer efficiencies of 85% are quite possible (7). The 

intrinsic coupling efficiency, n is the ratio of the light power accepted 
by the fiber to the power incidenton the fiber’s front facet. Thus n ^ 
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Figure 2. Cross-sectional views of the intensity distribution at the plane P. 2a is the intensity pattern for the rectanc^lar aperture; 
2b is for the circular aperture; 2c is the combined intensity pattern on the y,, axis; 2d is the combined pattern for yo= 
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Fiber 


Figure 3. Coupler efficiency defined as the product of the 
transfer and intrinsic coupling efficiencies. 
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depends on the accuracy of the hologram’s Imaging of the fiber's front 
facet. Leite and Soares (8) have shovm that values of n i above B5% can be 
achieved with only moderate difficulty. 

Taking these reported values for the respective efficiencies, it should 
be possible to attain coupling efficiencies, r\ of approximately 75%. 

IV. DOUBLE HOLOGRAM TECHNIQUE 

The single element holographic method for coupling, while Improving the 
coupling efficiency, still suffers, like the other coupling schemes, ^rom the 
severe alignment tolerances. Unlike the other coupling .jchemes, alignment 
tolerances can be relaxed by using a slightly revised holographic method. 

In this revised holographic method, the single hologram is split into 2 
holograms (Figure 4). The first hologram is designed to convert the non- 
symmetric diverging LD emission into a plane wave. The second hologram is 
designed to convert an Incoming plane wave into the fiber’s far-field 
diffraction pattern (Airy disk). By utilizing the fact that holograms have 
all of the pattern’s information stored in a small area, the alignment 
sensitivity, with respect to each other, of all four the components (LD, 
holograms 1 and 2, and the fiber) is greatly decreased when compared to the 
sensitivity of the other coupling methods.* 


*Lateral displacements between the coupler sections have virtually no 
consequence. However, tilt and rotation alignment become sensitive ((8) 
discusses methods for minimizing this sensitivity). Longitudinal displacement 
between the coupler and the source (LD or fiber) is believed to be no more 
sensitive than in the traditional techniques (20 microns for 1 dB loss) (4). 
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Figure 4. Double Holograun Configuration 


V. TWO-BEAM HOLOGRAPHIC FABRICATION 


Fabrication of the holographic couplers consists of merely making the 
hologram(s). For the single element device, the holographic recording schema 
employed is depicted in Figure 5a. The reconstruction (coupling) 
configuration is shown in Figure 5b. In the recording case, light emerges 
from the LD and Interferes on the photographic film/plate P with light from 
another laser focussed at the spot where the fiber's end facet will be 
placed.* The plate is then developed and reinserted into the system at plane 
P. When the hologram is illuminated by the LD radiation, the fiber’s 

coupling wave, U^, is generated which propagates Into the fiber. 

In the case of the double holographic element technique, two holograms 
are generated. The recording procedure consists of using plane reference 
waves, Upj and Up 2 » making complementary angles with the plane of the 
holograms in conjunction with the LD and fiber emissions (Figure 6). \7hen 
hologram HI is illuminated by the radiation from the LD, the plane wave is 
generated. This wave, in turn, generates the coupling wave when it 
illuminates hologram H2. 

VI. COMPUTER GENERATED HOLOGRAMS 

A. INTRODUCTION 

Another method of fabricating holograms is to use a computer and a 
plotter. In this method, the two-dimensional Intensity distribution resulting 


♦Notice, in Figure 5A, that the fiber is positioned conjugate to where it will 
be when LD light is coupled into the fiber (Figure 5b). 
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Figure 5b. Reconstruction (coupling) method for the single holograph! 
element coupler. 
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Figure 6. Recording scheme for the double hologram coupler. 
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from two beams interfering at some plane is calculated. This pattern is then 
converted into a form suitable for graphical representation by a discrete 
plotting device. A transparency Isfomed of the plot photographically or by 
some other method. The transparency then serves as the hologram. 

Calculation of the Intensity distribution resulting from the Interference 
of two beams Is a straight forward process. The diffraction patterns of the 
two radiation sources are known, or calculated, for the diffraction region In 
which they will be used (Fraunhofer or Fresnel). These fields are then 
propagated to the plane of Interference and summed.* The Intensity 
distribution of the Interference plane Is subsequently calculated. 

At this point the hologram's Intensity variation Is known by the 
computer. This pattern must then be displayed In a form such that a 
transparency may be formed. The simplest method of displaying the pattern is 
to plot It on a hard-copy graphics device. The resulting printout may then be 
photographed to form the transparency. 

The problem becomes one of how to accurately represent the Intensity 
pattern, consisting of phase and amplitude Information, by a discrete plotting 
device. 

The most widely used method for displaying phase and amplitude 
information on a discrete plotting device Is Lohmann's technique developed 
explicitly for computer generated holography (9). A binary (light-dark) 


*Thls plane defines the position between the two devices where the hologram 
will be placed for reconstruction. 
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plotting scheme is used since It eliminates the need for generation of gray 
levels. Another benefit of the binary scheme Is that In the photographic 
reduction step, development procedures are far less critical than for 
continuous-tone objects. 

Lohmann's method utilizes the fact that the phase of the offset reference 

wave, also used for reconstruction, varies linearly along the x axis.* 

Consider the system depicted In Figure 7 where a binary mask has been placed 

In the hologram plane. The mask Is primarily opaque with holes centered on a 

regular grid at the sample points. The linear distance Ax Is chosen such 

o 

that the phase of the offset reference wave Increases along x by 2tt radians 
through that distance. If two adjacent holes were spaced by Ax^ , then all of 
the light passed by the two holes will be of the same phase. If, however, one 
hole Is displaced a distance 1.25 Ax^ from Its adjacent one, as with holes C 
and B In Figure 7, then the light passed by that hole (C) will differ from the 
adjacent hole (B) by 2II+II/2 » 511/2 « II/2 radians. Destructive interference 
will occur between the light from these two holes. This technique of 
displacing the mask''s holes with respect to one another, in accordance with 
the spacing Ax^ , and thereby encoding the phase of each sample point, is what 
Lohmann devised. Displacement of the holes is only performed In one direction 
thereby Implying that therjs be constant phase along the orthogonal axis. The 
size of the hole corresponds to the amplitude of each sample. 

Choosing the linear displacement, Ax^ , which corresponds to a 2ir radian 
phase change of the reference wave establishes the recording, and thus the 

*This procedure is applicable for a planes reference wave. This technique 
would then be suitable for the two-element holographic coupler system, but not 
for the single element system. 
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reconstruction, angle for the hologram. Figure 8 depicts a plane wave 
incident cn a grid at an angle 0^^ . The path difference between the ad.lacent 
slits is given by: 

[AOPLJ - dsln©^ (13) 

The corresponding phase difference is: 

d -(•^)dsln0^ (14) 

We have chosen the linear distance ■ d to correspond to a radian phase 
^difference. So: 

2n - AX^slnG^ (15) 

The recording and reconstruction angle is then given by 5 

o 

The linear displacement distance, £^x^ , Is chosen based on, among other 
things, the resolution of the plotting device. Its maximum allowable plotting 
size, and the number of discrete points In the sample. 


B. FOURIER TRANSFORM HOLOGRAMS 

Lohmann first examined Fourlsr transform holograms (i.e. holograms of the 
Fourier transform of the object rather than of the object Itself) (9, 10). 
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The method used to generate the hologram relies on the plotting scheme 
described earlier. The procedure for producing a computer generated Fourier 
transform hologram la the following. The obiect, If It Is represented by a 
continuous function, Is decomposed Into discrete samples so that a fFT may be 
performed on It numerically. This sampling Interval should be chosen In 
accordance with the Whlttaker-Shannon sampling theorem thereby ensuring that 
the maximum amount of Information about the object Is retained. The Fourier 
transform of the object Is computed using a discrete FFT routine. Each sample 
of the Fourier transform of the object Is encoded via the Lotimann method. The 
plot Is photographically reduced, negated,* and reconstructed. Reconstruction 
Is accomplished by using a lens to perform an optical Fourier transform of the 
Image formed by the hologram (Figure 9). 

C. EXPERIMENT ~ A COMPUTER GENERATED BINARY FOURIER TRANSFORM HOLOGRAM 

A binary computer generated Fourier transform hologram of the two- 
dimensional character P was produced. A Digital Equipment Corporation VAX- 
11/780 computer with an associated Trllog printer/plotter was used to perform • 
the computations and plotting. The procedure followed In fabricating the 
hologram Is that outlined earlier. The character was defined In a discrete 
array of dimensions 64 x 64. By defining the object discretely, the need for 
accurate sampling of a continuous function was avoided. The Fourier transform 
of the object was performed by the IBM scientific support FFT subroutine HARM. 


*The aperture representing the phase and magnitude of each sample is darkened 
In. The use of the photographic negative as the hologram reverses the plotted 
contrast thereby making only the aperture transparent. 
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Reconstruction system for a Fourier transform hologram 


Each element In the Fourier transformed ob.iect has Its phase and 
magnitude encoded as an aperture within a square cell. The height of the 
aperture Is proportional to the magnitude of the sample, while the position of 
the aperture within the cell Is determined by the sample's phase. The size of 
the square encoding cell, along with the plotter's resolution, dictate the 
quantization levels for the phase and magnitude of the sample. Figure 10 
depicts two adjacent cells In which the phase and magnitude of two adjacent 
sample points of Fourier transform of the object will be plotted. The object 
was represented by a 64 x 64 element army. The Trllog plotter can plot 1320 
points across its 30.48 cm (12 Inch) plotting region. The plotter's 
resolution Is therefore 0.0231 cm. At this point, It would be possible to 
quantize the phase of each element Into 20 parts (1320/64 ■■ 20.62). While 
this would allow very good phase representation, the aperture would have a 
width of 0.0231 cm. This would result In a very dim Image because of the 
narrow aperture. The phase quantization level was chosen as 6 so that each 
phase region would be represented by 3 pixels.* The aperture width was thus 
0»0693 cm. The magnitude quantization level was net limited by the plotter, 
only by the length of plotting paper available. The magnitude quantization 
was chosen to be 7. Thus the height of the aperture could range from 0.0 to 
0.1386 cm. The linear extent of the square aperture cell (Figure 10) was 
0.1386 cm. Since the linear extent of the square cell corresponds to Ax^ , 
this Implies that Ax^ ■ 0.1386 cm. Ax^ corresponds to a 2tt radian phase 
change so by equation (16) we see that the recording and recontructlon angle, 
for He-Ne laser light (X - 6328 A) Is 0.0262 degrees. 


♦While this coarse phase quantization results in some degradation In the 
reconstructed Image, the image brightness was deemed more Important at this 
time. 
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The method for encodli|g the object’s Information Is to calculate the 
phase and magnitude of the Fourier transformed object and see what 
quantization bln each element's values fall Into. The phase quantization 
level of 6 means that the bins span the phases O-H/3, H/3-2ll/3, 2II/3-II 
II-4II/3, 4n/3“5ll/3 5n/3-2n, The object element being plotted has Its 
phase scanned and the appropriate phase bln noted. The maxlmun magnitude In 
the entire object array Is found. This point will have Its ent:l'j.‘e aperture 
filled In and Is therefore scaled to the magnitude quantization level of 7. 
The magnitude of each object Is then scaled to the maximum magnitude value. 
The magnitude of each element Is examined to see how much of the aperture 
should be filled In. The appropriate magnitude quantization amount Is 
noted. By preceding through the entire object array, the plot array Is 
determined and plotted. 


Figure 11 shows the original object plotted on the phase-magnitude 
grid. Figure 12 shows the Fourier transform of the object plotted on the 
phase-magnitude grid. Notice In Figure 12 that the sizes and locations of the 
apertures differ for different object points. 

Initial plotting of the Fourier transformed object showed that the FFT 
has too large of a dynamic range for quantization into only seven levels. 
Therefore, the logarithm of the FFTed object was used for magnitude 
quantization. To avoid arguments of zero, a one was added to each value In 
the FFT object. The side effect of this operation Is to enhance edge detail 
in the reconstructed image, as low frequencies are suppressed with respect to 
higher ones (11). 
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Original object plotted on the phase magnitude grid. 
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Fourier transform of object plotted on phase-magnitude grid. 


A listing of the computer program, CGHBLACKP FOR, used to generate the 


binary plots Is given in Figure 14. 


VII. CONCLUSIONS 


The use of holographic couplers for LD - fiber systems results In reduced 
alignment tolerances and higher efficiencies than the traditional techniques 
used In coupling. 

Fabrication of the couplers could be achieved via conventional 
holographic manufacturing techniques or by the computer generated approach. 
While both techniques have their advantages, overall the conventlally 
fabricated hologram Is the better of the two when used as coupling devices. 

The diffraction efficiency from a conventional hologram Is higher than from a 
binary computer generated hologram (12). Considerable higher resolution is 
obtainable In the conventional fabrication process than In the computer 
generated process. Holographic films have resolutions of 1250 llnes/mm (O.H 
microns) (13) while plotting devices resolution Is typically 200 microns. 

While binary computer generated holograms can be fabricated for non-exlstant 
devices, a mathematical description of the ob.lect Is necessary and may be hard 
to Implement on a computer. 

Each holographic coupler is fabricated for the particular LD - fiber pair 
with which it will be used. Spatial variations In the LD output beam change 
from LD to LD and therefore may not match the recorded LD beam pattern. 
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There are some limitations in the use of holographic coupling devices. 
Recording film sensitivity Is quite low in the 830 nm wavelength region. This 
effect requires relatively long exposure times over which the recording system 
must remain very stable. Fluctuations In the LT) output beam due to 
temperature Instability will Introduce a mismatching factor Into the coupling 
system. The hologram was developed for one particular LD beam pattern; when 
the reconstruction beam pattern differs from the recorded beam pattern, LH - 
fiber matching will degrade. 

Overall, the advantages of holographic couplers over the traditional 
couplers for LD - fiber systems warrant experimentation Into the hologram’s 
practical usefulness In coupling. 
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C Pro cram: CGHb lack? . FOR 

G 

G Tilts proi^rara gonerates a Fourier transform plot suitable for 

c conversion into a ho loiarram follovrlng pho tojfraphic reduction. k 

e The Trilo|r pr luter/plo t ter is used for generation of the binary 

. c plot, 

e 

e (GGIIblackP plots only black characters. the colored phase-innt;n i tiide 

e mrid Is not ar^uerated in this program). 


c author. P.Fuhr 

c date: March ~ April 1982 


complex a(64.64.1) 
complex b( 64,64) 

intei;er iplo tm( 64) , ipio tp( 64) 
dimension inv( 16) ,s( 16) 
dimension piz(6) 

integer rat 3) 
integer qz 

byte plot (222) 


common /tr idat/plo t , mco lor 


data ra/6,6,0/' 
data pl/3. 1415926/ 


larray containing complex image 
! * * magnitude and fjhase of 

ithis image's FFT 
fphase and magnitude plot arrays 
{arrays used by the FFT rout Ine IL4RT'I 

!PI array ranging from -(2/G):{:pi to 
!pi itself 

{another array used by HARM 


{array corresponding to the prograocned 
{plot data lines, there are 6 plot 
{lines for each of the 64 array rows 

{pass the data to the TRILOG plot 
{routines via a common block 

{size of the 3 dimensions of 'a' 

{PI 


c 


2 

ce 


3 


7 


£ 


t i 


init piz array 

inc>pi/3 . {PI increment (for piz) 

do 2 i? 1,6 

piz( i) * “((65. /64.) * pi) + (i^pizinc) 

before anything else, go grab the TRILOG 
call tri^init 

do 3 J> 1,220 {set bit 6 in each output data byte 

plot( J)» * 100 
cont inue 

plot(221)>5 {Hex 5 (Plot mode) 

piot(222)3l0 {Hex A (line feed) 

do 7 13 1,64 
do 7 J*l,64 

a( 1 , J , 1) 3c mp Ix( 0 . 0 , 0 . 0) {init the image array 

b( i • j) 3cmp lx( 0. 0, 0 . 0) {and the mag-phase array 

cont inue 

the character *P' 

define the character in the array 'a' 

do 20 131.3 

do 10 J3 1.20 

a( i.j, l)3cmplx( 1. 0,0.0) 
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le 

20 


25 

30 

40 


45 

60 


65 

c 

BO 


c o n t i iute 
cont time 

do 40 ls4. 10 
do 25 js|,5 

a( 1 , J • 1) -cmp Ix< 1 , Ot 0. 0) 
cont tnue 

do 30 js 16,20 

a( 1 , J , 1 ) -cinp lx( 1 . 0, 0. 0) 

cont inue 

cont Inue 


do 60 1=11, 14 
do 45 J=l,20 
a( l,J,l)-cmplx( 1.0, 0.0) 
cont inue 
cont inue 


do 80 t= 15.24 
do 65 J= i , 5 

a( i,J, I)=cmplx< 1.0, 0.0) 
cont luue 

cont inue 


Mhe charncter P have been defined in 'a,' 


c 

c 

c 

c 

c 

c 

c 

c 

c71 

c 

c 

c , 
c 

c73 

c 

c 

c74 

cc 

c 

c75 

c 

c77 

c 

c 

c 

c 

c 


88 


GGllOLAGK producea output in btack only 
tncolor^i fyellow ribbon 

call trl—color Ichangre ribbon 

no need for the grrid 

da 77 k= 1 ,64 
do 71 1=1. 195 
pIot( l)= "177 

call tr i-dump 

do 73 i= t .220 
piot( l)= "100 

do 74 1=2, 195,3 
plot< 1)= "120 

do 75 1= 1 . 1 1 
call tri_durop 

cont inue 

mco lor=3 
coll tr i_co lor 

ca 1 I tri — top 

t>T»e *, • enter o to plot origrinal. f for FFT, i for Fiiiv(F)’ 
accept 80. qz 
forma t ( la 1 ) 

lerr = 0 Jinlt HARM error flag’ 


Tevory bit on 
! p 1 o t this 1 1 lie 

! every bit off 

levery 6th horizontal space is a vertical line 
Ki.e. every 3rd byte ="120) 

111 times 

levery array row 

fblue ribbon 
Ichange ribbon color 
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e 


83 

e 


c 

as 


c 

c 

e 

e 

90 

95 

°00 


1005 
107 
1 10 


lf:set*1 fhov» HAHj'I greiicrnte 1 1» own interiiMl »in« 

land co»l|ie tirbl## for the FFT 

1 f ( qz> eq« *0 ' >iro to 85 flf requested, don’t do the FFT 

call HARTK a»m, inv, a . I fae t « ierr) !SSP FFT routine 

If ( ierr .eq.Ol vo to 83 

t>T>e ’ FFT error code »*.lerr 

1 f ( qz . ne • ‘ 1 ’ ) KO to as Idon* t perform Flnv uiileaa reqiiealed 

ifaet*-l Mtave HAIU1 greucratc Ita own internal atna 

land coainti tab lea for F liiverae 
call 'flAlUK a , B, Inv.B . I fae t , lerr) !SSP FFT routine 
lf( ierr<eq.O)$roto 85 

type (•,’ inverae FFT error code *',lerr 

determine the magnitude and phnae of each object element 

do 100 1*1.64 
do 100 J*1.64 

amag*caba( a< 1 , J . 1 ) ) Ifind the element’a magnitude 

If ( rea 1 ( a( 1 , J , 1 > ) . ne . 0. 0) iro to 90 Idoes the real part equal 0? 

x»pi/'2. 

y*a Ireaffl a( 1 • J . 1 > ) Ij^et the iinag'tnary componeat 

apha*a lfn( s. y) 'the phnae * pi/2 I imea the a iffn (+ or **) of 

1 the hnagrlitary part. (If the real component of 
Ithia iiuuLer >0, then the phase anirlc 1* * or 
1~ pi/'2. if the Iioaeinary part la uegative 
1 then (he phase an^le is negrative. 

firoto 95 


aplia*a t un2( a iiuagrt a( 1 I J . 1 ) ) . rea 1 ( a( i I J . 1) > ) luel^yrraine the pliaae anffle 


b( 1 , j) *cmplx(amasr,apha) 
cont inue 


Iput the mag;nitude and phase into the real 
land iina grinary parts of b 


type *. • do you want to s tore the mag-phase data (y or n)?' 
accept 88. q7 
1 f ( qz . ne . ' y’ ) go to 104 
do 110 i«l,64 
do 107 j * I > 64 

zmag*rea 1 ( b( 1 . J > ) Ireal part ( r.iasrn 1 1 ude) 

zpha*a Imag(b( i , j) 1 Icompiex part (phaseJ 

wr 1 te( 9 . 1005) 1 , j . v.inag, zpha 

formate bC . 12, 52, ’)» * ,e I5.8.0x,e 15.8) 

cont inue 
cont inue 


c 

104 

c 


» 

t 


120 


find the maxiffliim magnitude value in the b ar.-ay 

zmax*log(real(b( 1 , 1) )^1 .0) Intake the 1st element t!ie largest (start) 

take logarithm because FFT has too large of a dynamic range 
do 120 1*1.64 
do 120 j*1.64 
z*log(real(bt i,j))vl.0) 
zmax*aiaaxl ( z.zmax) 


cont! nue 


Icheck every element 

luse the maximum value function for the 
Icompar ison 
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e 


e 


e 




c 


2oe 

cc 

cc 

c 

c 


c 


c 


400 

c 


405 


there are 7 levela of maifnitude quantization 
zinIevl*;imaJc/7.0 

there* are 6 levels of quantization for the phase 
a£plevl*t> 

determine ihe phase quantization of each element 
do 500 ;«1,64 


do 200 J>1,64 

p*a inui|f( b( J ) ) !get the phase 

if(p, ll.pIzC 1)> lplotp(j)>l 
lf(p.gc.plz(l).and.p.lt.piz<2>)iplotp(j)>2 
lf(p.ge.ptz(2) .an'd.p. lt.piz(3>) iplotp(j)*3 
lf(p.ge.piz(S).aud.p. lt.piz<4>) iplotp(j)*4 
if(p.ge.plz(4).and.p. lt.piz(5)) iplotp(j>*5 
if(p,ge.plz(3)) iplotp(j)e6 


!de termine 

! which 

fsegtnent the pliase 
tshould be encoded in. 


determine the proper amplitude quantization level for each element 
amag* log( rea 1( b< 1 1 J ) ) + 1 . 0) fget the magnitude 

z leve I'amag/zmlevl Iscale the data 

Iflzlevel. It. 1.0) lploira(j)s0 

If (zleve 1 . ge . 1 . 0. and . z le ve 1. It. 2.0) Iplotmtj)*! 
if (z levs 1 .ge .2.0. and .z ieve 1. It. 3.0) ip]otm{j)32 
if (z levs 1 . ge .3. 0.aitd .zlevel. It. 4.0) iplotm(j)=3 
if ( zleve 1 . ge .4.0. and .zleve 1. It. 5.0) lplotm(j)>4 
1 f (zleve 1 .ge .5.0. and .z leve l.lt.6.0)iplo tm( j) *5 
If ( zle ve 1 . ge . 6 . 0) ip lo tm( j ) «6 !the largest magnitude 

cont Inue 


nez-t. map this plotting data into the output plot Hues (6 plot 
lines per 1 array line). 


do 410 k>l,6 
lplevl=7-k 
do 400 j-1.64 


scan a plot 1 ine 

1 f ( Ipio tm( j ) . 1 1 . Ip le vl ) go to 400 Igo 

there is something to be plotted on 
l = 3*j 

lf( iplotp(j) .eq. l)plot( 1+I)s-170 
lf( iplotp( j) .eq.2)plot( 1+I)» *107 
lf( lplotp( J) .eq.3)plot( l-4-2)» *170 
if ( iplo tp( j) .eq.4)plot( 1+2) = *107 
if( lplotp( j) .eq.5)plot( I+3)«*170 
if( lplotp( J) .eq.6)plot( 1+3) » “107 

cont inue 


on to next array element 

this line, find the phase position 


plot out this line 
call trl...dump 
call tri_dunp 


!2 vertical lines per amplitude quantization 
! level. (12 Hues per array row). 


do 405 kl> 1,220 

p lo t ( kl ) = * lOO treset the plot array 

cont Inue 
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410 continue 

500 continue 

e top 
end 
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C ?nOGF’fih NAHE5 FF.FOR 

c 

c This p-roarani calculites Ih# intansily for two electric fields 

c inlerfarina it s ^lan* i distinct z between them* Fields are 

c Cslcuizted for a rectangular srerture and a circular s*'erture. 

c T.-iSs* fields art then coRibined at the holosrafhic (interference) 

c ^lane» The resultins intensity psttern is subseauentN scanned 

c for its {iiaNijftuni value* Trte output data is scaled such that the 

c ir<s>^iniuni value will Just ssturate the output plot* F'lottins is 

c performed on a DEC VT-11 display and a t'ersatec printer/plotter . 

c 

c . author t P. Fuhr 
c date: April 1982 


FEALU AVAL< I024>rout>{( l024)fOut«< 1024) 

DIMENSION IBUF(aOOO) 

L0GICAL4C1 STRNG(21)rae 
rtall(4 Ixf IVf lUembdafSrpifKz 

CALL INIT( IBUFfoOOO ) linit graphics diss-lay 

c ecGuire user input 

type enter xsero endpoint (in microns)' 
accept iTrxoend 

type tr' e symmetric sraph will be generated' 
xostrt * l.e-6 {start near zero 

c we will senerate 512 points the other 512 points come from the 

c symmetry of the system 

«Dinc * (Koend - >:ostrt )/512. 

ty«-e #>' enter constants (all lengths are in microns)' 
type ¥»' enter x dimension of rectangular aperture' 
accept 

type I*’' enter y dimension of rectangular aperture' 
accept ifly 

type *r' enter diameter of fiber core' 
accept *rl 

type >».f' enter wavelength of light' 
accept If lambda 

type If' enter distance betwenn apertures and interference plane 
accept IfZ 

c all lengths are in MICRONS 

pi*3.141SS 

Kz=( 2*|pi )/lambda !propssstion constant 

yzero * (pi * iaiTibds I z) / (2* * ly) 
xzsro=>;oitrt 

noff=5l2 ! internal array offset 
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TYPE tf' in?ut Ysaro vslut (y or n)' 
scce^'t 987fao 

if( Qa.n».'y' .jm*i*c»a»ne»' Y' )^oLo 87 
type tr* enter Ysero' 
sccept .t>ysero 

87 do 600 isl»512 kalculate the intensity at each point 

type Jitfi !let the* operator Know whst point you're on* 

call rect( ansuer» 2 f l5-{f lyf laiTibdafXserofyzero ) !rect aperture' 
call circ< answr2»Zf IfKsfXserofysero) Icircular aperture 

first » ansuer I ansuer isouare the field 

second * answr2 $ ansur2 Isauare the field 

third » 2* % ansuer i ensur2 ’cross ten 
8val( iriioff ) » first 4- second f third Isun ail components 
>:sero ■ xzero + xoinc ! increment output field position 

600 continue 

c now fill the bsainnina of the array 

nr, *1024 

do 605 i*l*S12 
aval( i >*aval( nn-i ) 

60S continue. 

N3*0 

C FIND HAX. VALUES IN ARRAYS 

AHAX-0 . 
smin ■ l*elO 
DO I K*lrl024 

IF ( ABS(AVAL<K)).GT,Af1AX) A«AX*ABS( AVAL( K ) ) 
if( 3val( K )*lt»sniin ) smin = avsKK) 

1 CONTINUE 

type print data (y or 
accept ?37fQa 
9S7 formaU lal ) 

if( ao*ne*'y' *^nd*OQ*ne*' Y' )soto 16 
do 73 i*lfl024 
print 933f i rsvaK i ) 

933 format< 2xf i4r2x»ei6.8 ) 

73 continue 

print 934»ama}Cf3min 

934 formsU' AhAX* '*el6.3f' AMIN* 'rel6*8) 

C SIZE THE SCREEN f DRAW BORDER 

16 .CALL SIZE(7*5»7.S) 

CALL SCAL(0.f0.»72.f72. ) 

CALL AXES(0.»0.»72.f72. ) 

CALL PLOTl 

if( smsK.ne^O* )5oto 17 
amax* 1 » 
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c i:ale Iht data to AHAX and AJ'.IN 

17 'sscile » 70./81.1X 

18 ^scaUsS0./1024. 
do 23 i*ift024 

outy( i )s3vsl( X )'’vsc3lf I'XinH forsst th? w off sat 

outx( i )af lo 2 t( i )t>i5cal« flO. , !o(^ tne ),* a rxa offs-et 
23 continut 

c flot out thf data 

DO 30 i»lrl024 
K*out«< i ) 
y»out«( i ) 

ciU 8?'nt(M)-yr0f8) 

30 CONTINUE 

C LABEL GRAPH 

TYPE 3 

3 FORHAK ' X-AXIS LABEL?' ) 

CALL GET37R- S»STRNCf20fERR) 

#1*0 I horizontal Isbai 

X»40. 

Y*7, 

CALL FLDT2(STRNGfXYY.!ti) 

TtPE 4* 

4 FORMAT^ ' Y-AXIS LABEL?' ) 

CALL GETSTR(S»STRNG*20fERR) 

X*1C . 

ii*l Ivtrtical isbai 

Y^'SO^ 

CALL PL0T2' STRNGfXfYfn.) 

TYPE S 

5 F0RHAT(Xf'C0PY; O^sNOr 1=YES?' ) 

ACCEPT *f N1 

iF<Ni noas.u 
11 CALL PLOT 

CALL DONE 

15 STOP 

END 


SUBROUTINE TQ DRAW BOARDER 


SUBROUTINE PLOTl 
call a^n t( 1 0 » f 0 • » 0 > — 0 ) 
call vect( 0*f70. ) 
call sf^nt( 70i rOf fOf-8 ) 
call vacl( 0. f70. ) 


.•lower left 
!left vertical 
•lower risht 
Iriaht vertical 
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call BF-nU 10* ?0, >0>-3 ) 
call vecK 60. *0* ) 
call B.*nU 10. y.’O* fOr-B ) 
call vect( 60. fO. ) 
call 10. »35. f0»-8 ) 
call vecU 60. rO. > 
call 3.'='nl( 40. fO. >0»-3 ) 
call veci( 0. f70. ) 
return 
END 


.'lower left 
ibottok 
•upper l«fi 
! top 

jriidcfle left 
Ihorizontal axis 
liTilddle bottoe 
.'veritcal sxis 
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SUl '•routine to draw text 

SUBROUTINE PL0T2( STRNG»X»Y.r» ) 
LOGICAL*! STRNG<21) 

L*LEN( STRNG ) 

IF(M.EQ.l ) GO TO 22 
HORIZ. TEXT 
CsX-L/2. 

CALL HTXT(C»Y»LfSTRNG) 

GO TO 24 
UERT. TEXT 
C=Y-L/2 

CALL APNT(X?C»0^-B) 

CALL VTXT(XrCfLfSTRNG) 

CONTINUE 

RETURN 

END 


subroutine rect( answer fZflxflyflatribdsrxzepofyzero ) 
iniplicit real*4 (s-ifO-z) 

compute the far-field due to a rectangular aperture 

calculate the amplitude constants 

if( Ix.ne.O.O.and.ls.ne.O.O )Soto 5 

answer =0.0 

return 

a « (lx » ly) / (lambda * z) 

if ( yzer 0 . ee . 0 . 0 )yzer o=l . e-5 

sincy = sin(yiaro) / yzero 

xOara s (lx * xzero) / (lambda * z> 

if( xOars.ee.O.O )x0ars=l .e-5 

sincx = sin(x0ar2) / xOsrs 

answer = a * sincx * sincy 

return 

end 

subroutine circC answr2>z?l»Kz»xzerofyzero ) 
compute the far-field out of a circular aperture 
implicit reai*4 (s-z) 
if( l.ne.O.Ovand.Kz.ns.O.O )soto 5 


CJ CJ U t_> CJ O O O O '-J CJ ' J o C.» iJ O CJ C J t J 


c 

c 


1 


3 

4 

10 
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b-<Ks*ll:l)/(8.»2) 

K12z s <Kz * 1) / (2. * 2 ) 

tessra a (xzero * >izero) + (yzsro » yzero) 

besarg » K12z t sarKbesarg) 

angur2 * (2* 4 b / bea-arg) * bssseK beaarg ) 

return 

end 


function bssael(x) 

compute the Bessel function of the first Kind - first order 
b» using the scientific subroutine pecKsged program B£S-J*F0R 


z*>{ 

.■argument 

of the Bessel function 

n=l 

list order 

Bessel function 

da.Ol 

! accuracy 

desired 


call besJ( ZfHfxJ^df ier ) 

if( ier.ea.OSoto 10 IchecK for errors 

bessel * 0. 

goto( 1 f 2f 3f 4 )ier Itype the corresponding error message 

type *»' order of Bessel function less than zero' 

return 

type *f' Bessel argument negative or zero' 
return 

type >f' reauired sccuracy of 0f007 not obtained' 
return 

type *»' range of Bessel order not correct < see S3P doc p*3-35)' 
return 

bessel * >?J Ireturn with the answer 

return 

en*:l 


SUBROUTINE BESJ 


PURPOSE 

CDHPUTE THE J BESSEL FUNCTION FOR A GIVEN ARGUhENT AND ORDER 
USAGE 

CALL BESJ(X»N?BJfDfIER) ! 

DESCRIPTION OF PARAMETERS . 1 

X -THE ARGUMENT OF THE J BESSEL FUNCTION DESIRED I 

N -THE ORDER OF THE J BESSEL FUNCTION DESIRED I 

BJ -THE RESULTANT J BESSEL FUNCTION . | 

D -REGUIRED ACCURACY I 

lER-RESULTANT ERROR CODE WHERE [ 

IER=0 NO ERROR I 

IER=1 N IS NEGATIVE I 

I 

• i 

ii 
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IcR»2 X IS NEGATIVE OR ZERO 

IER»3 REQUIREH ACCURACY NOT 0&TAINEH 

IER»4 RANGE OF N COMPARED TO X NOT CORRECT (SEE REMARKS) 

REMARKS 

N MUST BE GREATER THAN OR EQUAL TO ZEROr BUT IT MUST BE 
LESS THAN 

20fl0»X-X*» 2/3 FOR X LESS THAN OR EQUAL TO 15 
904X/2 FOR X GREATER THAN 15 

SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 
NONE 

METHOD 

RECURRENCE RELATION TECHNIQUE DESCRIBED BY H* GOLDSTEIN AND 
R.M»' THALER RECURRENCE TECHNIQUES FOR THE CALCULATION OF 
BESSEL FUNCT10NSSM»T.A*C*fV.||»PPa02-J08 AND I,A» STEGUN 
AND M. ABRAMOUITZf' GENERATION OF BESSEL FUNCTIONS ON HIGH 
SPEED COMPUTERS' »M.T*A,C.fV.Hrl957»PP. 255-257 


4 * 


4 4 > 4 « 4 » » 4 • 4 


SUBROUTINE BESJ( XfNi BJf Df IcR ) 
BJ**0 

IF(N)i0f20»20 
10 IER»1 
RETURN 

20 IF(X)30f30f31 

30 IER-2 
RETURN 

31 IF(X-15. )32»32r34 

32 NTES'T»20» + 10»4X-X*¥ 2/3 
GO TO 26 

34 NTEST»90.+X/2» 

36 IF(N-NTEST)40»33f38 

35 IER*4 
RETURN 

40 IER»0 
N1=N+1 
BPREV*.0 

COMPUTE STARTING VALUE OF M 

IF(X-5. )50f60f60 
50 MA»X+6* 

GO TO 70 

60 MA*U4*Xf60./X 
70 ME*N+IFI>!(X)/4+2 
MZERQ=MAX0( MAfhB) 

SET UPPER LIMIT OF M 
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HHAX=NTEST 

100 DO 190 H«hZER0f«hAXf3 
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SET F(h)fF(M-l) 

FH1*1»0E-2B 

FMa.O 

ALPHA=.0 

IF<H-< H/2)»2)120.110»120 
110 JT»-1 

GO TO 130 
120 JT*1 
130 H2=H-2 

DO UO K«1»H2 
MK»H-K 

BliK»2.»FL0AT(MK)»Fl11/X-FM 

FH*FH1 

FM1»BMK 

IF< KK-N-l)lS0»140fl50 
140 BJ=BMK 
150 JT»-JT 
S-l+JT 

160 ALPHAsALPHA+BMK*S 
BMK*2.*FM1/X-FM 
IF(N)ia0-rl70fl80 
170 BJ=BMK 
180 ALPHA=ALPHA+BHK 
BJ=BJ/ALPHA 

IF( ABS( BJ-BPREU )-ASS D4BJ ) )200r 200r 190 
190 BPREV*BJ 
IER-3 
200 RETURN 
END • 
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1* C*C* Tinflieriianr Opl*r 15» 2432('l?76)t 

2* H* Saruuaiari and K. NawataK» A^pI* Opl*r 1B» 1847(197? )» 

3* L*C* Cohen and H*V* Schneider? AppI* Opi«? 13? 1021(1974)* 

4* M* Saruwaiari^ and T* Susie? IEEE JGE? QE-17? 1021(1981)* 

5* private conversation (General Optronics Corp») 

it A* Jordan? 'Electroaasnetic Uaves and Radiatins S»stens” 

7* H*P» Jordan? L* Solumar? and P»St*J* Russel? lEE J* Hie* Opt* 
Acoust*? 2? 156(1978)* 

8* A»H*P*P Leitc and 0*D*D* Soares? Proc* S*P»Z*E*? 213? 10(1981)* 

9* B*R* Brown and A*U* Lohaann? AppI* Opt*? 5? 967(1966)* 

10* A*U* Lohmann and B*P* Paris? AppI* Opt*? 6? 1739(1967)* 
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